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Abstract. We study the possibility of finite-time blow-up for a two dimensional Broadwell model. 
In a set of rescaled variables, we prove that no self-similar blow-up solution exists, and derive some 
a priori bounds on the blow-up rate. In the final section, a possible blow-up scenario is discussed. 



1 - Introduction 

Consider the simplified model of a gas whose particles can have only finitely many speeds, say 
Ci,. . . ,cjv £ JR n • Call m = Ui(t,x) the density of particles with speed q. The evolution of these 
densities can then be described by a semilinear system of the form 



Here the coefficient a^fc measures the rate at which new z-particles are created, as a result of 
collisions between j- and /c-particles. In a realistic model, these coefficients must satisfy a set of 
identities, accounting for the conservation of mass, momentum and energy. 
Given a continuous, bounded initial data 





(1.1) 



Ui(0,x) = Ui(x) 



(1.2) 



on a small time interval t £ [0, T] a solution of the Cauchy problem can be constructed by the 
method of characteristics. Indeed, since the system is semilinear, this solution is obtained as the 



fixed point of the integral transformation 



Ui(t,x) = Ui(x - Cit) + y2a ijk UjU k (s, x - Ci(t - s)) ds . (1.3) 

Jo j,h 

For sufficiently small time intervals, the existence of a unique fixed point follows from the contrac- 
tion mapping principle, without any assumption on the constants a^-fc. 

If the initial data is suitably small, the solution remains uniformly bounded for all times [3]. 
For large initial data, on the other hand, the global existence and stability of solutions is known 
only in the one-dimensional case [2, 6, 10]. Since the right hand side has quadratic growth, it might 
happen that the solution blows up in finite time. Examples where the L°° norm of the solution 
becomes arbitrarily large as t — ► oo are easy to construct [7]. In the present paper we focus on 
the two-dimensional Broadwell model and examine the possibility that blow-up actually occurs in 
finite time. 

Since the equations (1.1) admit a natural symmetry group, one can perform an asymptotic 
rescaling of variables and ask whether there is a blow-up solution which, in the rescaled variables, 
converges to a steady state. This technique has been widely used to study blow-up singularities 
of reaction-diffusion equations with superlinear forcing terms [4, 5]. See also [9] for an example 
of self-similar blow-up for hyperbolic conservation laws. Our results show, however, that for the 
two-dimensional Broadwell model no such self-similar blow-up solution exists. 

If blow-up occurs at a time T, our results imply that for times t — > T— one has 

, , 1 ln|ln(T-i)| , 

IK*)|U > 5 l T -t • (L4) 

This means that the blow-up rate must be different from the natural growth rate ||^(i)|| L oo = 
0(1) • (T - t)~ l which would be obtained in case of a quadratic equation u = Cu . 

In the final section of this paper we discuss a possible scenario for blow-up. The analysis 
highlights how carefully chosen should be the initial data, if blow-up is ever to happen. This 
suggests that finite time blow-up is a highly non-generic phenomenon, something one would not 
expect to encounter in numerical simulations. 



2 - Coordinate rescaling 

In the following, we say that P* = (t*,x*) is a blow-up point if 

limsup Ui(t,x) = oo 

X — >X* , t—rt* — 

for some i <G {1, . . . , N}. Define the constant 

C = max \ci \ . 

i 

We say that (t*,x*) is a primary blow-up point if it is a blow-up point and the backward cone 

T={(t,x); \x-x*\<2C(t* -t)} 
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does not contain any other blow-up point. 

Lemma 1. Let u = u(t,x) be a solution of the Cauchy problem (1.1)- (1.2) with continuous initial 
data. If no primary blow-up point exist, then u is continuous on the whole domain [0, oo[ xM n . 

Proof. If u is not continuous, it must be unbounded in the neighborhood of some point. Hence 
some blow-up point exists. Call B the set of such blow-up points. Define the function 

<p(x)= inf {t + C\x-£\} . 

By Ekeland's variational principle (see [1], p. 254), there exists a point x* such that 

C 

<p(x) > tp(x*) — — \x — x*\ 

for all x € M 2 . Then P* = (tp(x*), x*) is a primary blow-up point. □ 
Let now (t*,x*) be a primary blow-up point. One way to study the local asymptotic behavior 
of u is to rewrite the system in terms of the rescaled variables Wi = «;» (r, 97) , defined by 

r = -ln(t* -t), 

v = e r x = ^JZ^l (2.1) 
1 t* -t 

Wi = e~ T Ui = (t*—t)ui. 
The corresponding system of evolution equations is 

d T Wi + (a + n) ■ V v Wi = -Wi + ^aijkWjWk i = l,...,n. (2.2) 

Any nontrivial stationary or periodic solution w of (2.2) would yield a solution u of (1.1) which 
blows up at (t*,x*). On the other hand, the non-existence of such solutions for (2.2) would suggest 
that finite time blow-up for (1.1) is unlikely. 



3 - The two-dimensional Broadwell model 

Consider a system on M 2 consisting of 4 types particles (fig. 1), with speeds 

C 1 = (M), C 2 = (l,-1), C 3 = (-l,-l), C 4 = (-l,l). 




The evolution equations are 

(d t ui + Ci • Vui = u 2 u 4 - U1U3 , 
t u 3 + c 3 • Vu 3 = u 2 u 4 - mu 3 , 
d t u 2 + c 2 • Vu 2 = «iu 3 - ii 2 u 4 , 

k d t U4 + C4 • VU4 = U\U 3 — U2U4 . 

After renaming variables, the corresponding rescaled system (2.2) takes the form 

' d t wi + (x + 1)9^1 + (y + l)d y w 1 = w 2 W4 - toiw 3 - w\ , 
d t w 3 + (x - l)d x w 3 + (y - l)d y w 3 = w 2 w 4: - wiw 3 - w 3 , 

d t W2 + {X + 1)9x^2 + (y - l)d y W 2 = WlW 3 - W 2 W4 - w 2 , 

. <9 t -W4 + (x - l)d x Wi + (y + l)d y W4 = wiw 3 - w 2 Wi - W4 . 



(3.1) 



(3-2) 



Our first result rules out the possibility of asymptotically self-similar blow-up solutions. A 
sharper estimate will be proved later. 

Theorem 1. The system (3.2) admits no nontrivial positive bounded solution which is constant 
or periodic in time. 



Proof. Assume 



< Wi(t, x,y) < k 



(3.3) 



for all t, x,y, i = 1, 2. Choose e = e /2, so that 



1 

£<-, 
K 



ee lKX <- a; €[-1,1]. 



Define 



Qu(t,v) = J [( 



(1 _ ee 2KX ) Wl (t, x, y) + (1 - ee- 2KX )w 4 (t, x, y) 



dx . 



Qi 4 {t) = sup Qu(t,y) , 
l»l<i 

Restricted to any horizontal moving line y = y(t) such that y = y + 1 (fig. 2), the equations 
(3.2) become 

d t wi + (x + l)d x w\ = w 2 W4 - wi^3 - wi , 

9tU)4 + (x — l)d x W4 = WiW 3 — W 2 U)4 — W4 . 

A direct computation now yields 
j t Qi4{t,y(t)) 

< -2ek 



< —£K 



e 2KX (l + x)w 1 + e~ 2KX (l -x)w 4 dx + {ee 2KX - ee~ 2KX )(wiw 3 - w 2 W4) dx 
1 L - 1 J-i 



1 r pO 



e 2KX {\ + x)w! + e~ 2KX {l - x)w4\ dx - en J e~ 2KX (l - x)w4 . 



+ x)wi dx+ I ene~ 2KX W4 dx + I eKe 2KX wi dx 



~ eK j e 2KX (1 + x)w! dx + J , 



< —EK 



e 2KX (l + x)w 1 + e- 2KX (l - x)w4 



dx . 




Call 



I(t,y) = / w 1 (t,x,y)dx. 



The definition of e and the bound (3.3) on w\ imply 

J e 2KX (I + x)w± dx > 2e j (l + x)w 1 dx 

r-l+I/n 

> 2e J (1 + x)k 



dx 



el 2 Ik 



From this, and a similar estimate for W4, we obtain 

J [e 2KX (l + x)w 1 +e- 2KX (l-x)w i jdx>^J Wl dx^j + w^dx^j >-J^- 

Since e < this yields 



j t Qi4t,y(t)) < - £ - Q 2 1A (t,y(t)) . 



(3.4) 



Observing that the Cauchy problem 



■ e 2 2 



z(0) = 4k 



has the solution 



by a comparison argument from (3.4) we deduce 

Qu(*)<(^ + ^ 
Since 

/ wi(t,x,y)dxdy < 4Q 14 (t), 
-lJ-i 

and since a similar estimate can be performed for all components to,, we conclude 

j J Wi(t,x,y)dxdy ^ 4 (^ + V"*) " (3 ' 5) 

The right hand side of (3.5) approaches zero as i — ► oo. Therefore, nontrivial constant or time- 
periodic L°° solutions of (3.2) cannot exist. □ 



4 - Refined blow-up estimates 

If (t*,x*) is a blow-up point, our analysis has shown that in the rescaled coordinates r, £ the 
corresponding functions must become unbounded as r — ► oo. In this section we refine the 
previous result, establishing a lower bound for the rate at which such explosion takes place. 

Theorem 2. Let u be a continuous solution of the Broadwell system (2.1). Fix any point (t*,x*) 
and consider the corresponding rescaled variables r, If 

max Wj(T,£i,£ 2 ) < #hiT i = 1,2,3,4, (4.1) 

l£i|>l£2|<i 

for some 6 < 1/4 and all r sufficiently large, then 

lim^(r,O = t = l,2,3,4, (4.2) 

T-»(X) 

uniformly for £ G iR 2 in compact sets. Therefore (t*,x*) is not a blow up point. 

Since to* = (£* — and r = | ln(i* — i)|, the above implies 
Corollary 1. If (t*,x*) is a primary blow-up point, then 

I / m **-* 1 
hmsup \u(t,x) ■ ; -t > -. 

s ^ x ., tit*- 1 V " ln|ln(t* -t)| " 4 

Proof of Theorem 2. 
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Let Wi = Wi(t,x,y) provide a solution to the system (3.2), with 

< Wi (t, x, y) < 6 In t = k(t) (4.3) 

for all t > to and x, y € [—1,1]. The proof will be given in two steps. First we show that the L 1 
norm of the components Wi approaches zero as t — » oo. Then we refine the estimates, and prove 
that also the L°° norm asymptotically vanishes. 

STEP 1: Integral estimates. Consider the function 



Quit, y) 



e 



2fc(t)(x-l)\ / e -2fc(t)(x+l) 



1 ) W!it,x,y) + ( 1 ) w±{t,x,y) 



dx 



with kit) as in (4.3). As in the proof of Theorem 1, let t <— ► yit) be a solution to y = y + 1. Then 
J t Q^ Vit)) =f[-i x ~ l)k'e 2k{t ^-^ Wl + ix + l)k'e- 2k ^ + V Wi \ dx 

rl / e 2fc(t)(x-l)\ 

+ / ( 1 j [-(x + l)w lx + w 2 Wi - W!W 3 - wi] dx 

rl / g-2fc(t)(x+l)\ 
+ J I 1 J [ — (x — l)w ix + W1W3 - W2W4 - iu 4 ] dx 

To estimate the right hand side, we notice that 

rl / 2fe(t)(x-l)\ /-l 

A = J ( 1 J [(1 + x)toi x + wi] dx > kit) J ix + l)e 2fc(x_1) ti;i 

/•l / p -2fe(t)(a;+l)\ rl 

B = J M J [(x- l)w 4 * + ^ 4 ]dx > fc(t)y (l-x^-^+^dx 

/•l / e 2k[x-l) e -2fc(x+l)\ 

c = y (^i^3 - ^2^4) ( — 2 2 — ) dx 

rl e 2k(x-l) rO e ~2k(x+l) 

< kit) / W\dx + kit) / w^dx . 

Jo ^ J_ x 2 

Therefore, 

^-Qi 4 (t, yit)) = t [ - (x - ljfc'e 2 *^^- 1 )^! + (x + l)k'e- 2k ^ x+1 ^w 4 ] dx-A~B + C 
dt J _i 

<k\t) J' [il-x)e 2k ^ x -^w 1 + il + x)e- 2k ^ x+1 ^w 4 ]dx 

_M f 1 [(l + x)e 2fc «( a; - 1 ^ 1 + (l-x)e- 2fc «( a;+1 V]dx. 
If & it) > 1/2, we claim that the following two inequalities hold: 

(1 _ a .) e 2fc(t)(x-l) < ! _ e 2fc(t)(x-l) ) 

(1 + x)e" 2fc(t)(a:+1) < 1 - e - 2fc (*)( a: + 1 ) . 



(4.4) 



To prove the first inequality we need to show that 

h k (s) = 1 - e 2ks + se 2ks > for all s € [-2, 0] . 

This is clear because hk(0) = and 

h' k (s) = e 2ks (l - 2k + 2ks) < s€[-2,0] 

if k > 1/2. Hence hk(s) is positive for s £ [—2,0], as claimed. The second inequality in (4.4) is 
proved similarly. 

When t > t = e 1/(26) one has k(t) > \ and hence 



dt Qu{t, y(t)) < k'(t)Q 14 - ^ J 1 ^l + x)e 2k ^-^w 1 + (1 -x)e- 2k ^^w,}dx. 



j t Qi4{t,y(t))<k'(t)Q 14 k{t) f 
Setting / = J_ x w\dx, we obtain 

r l ,-l + I/k(t) j 2 

/ (l + x)e 2k ^- 1 ^w 1 dx > / (l + x)e- ik ^k(t)dx = e" 4 ^ -f— 

J -i J -i 2k(t) 

Using the above, and a similar estimate for the integral of w A , we obtain 

^ C [(1 + x)e 2k ^ x -^ Wl + (1 - x)e- 2k ^ x+ ^w±]dx 
2 J-i 



> 



e -4fe(t) 

"4 

e -4fe(t) 



u>i(ix + I / u>4<ix 



(4.5) 



> — — Q?4 



Calling 

from (4.5) we deduce 



Qu{t) = max Qu{t,y) , 



J -4fe(t) 

^0i4(*) < k'(t)Q 14 (t) — Q 14 (t) 2 . 



Recalling that k(t) = 8lnt for some < 6 < 1/4, the previous differential inequality can be written 

as 

f t Qu < \q» - ^Q? 4 . (4.6) 

Notice that Q\4(to, 2/(*o)) < 2fc(t ), and define the constant 

A = max{2£;(to)to" 4e > 8 ( x - 3 #)} • 

Then the function 

z(t) = A t ie - 1 

satisfies 

j t z{t) >-z- -^z 2 z(t ) > Qu(to,v(to)) • (4.7) 
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Comparing (4.6) with (4.7) we conclude 

Q 14 (t)<z(t) t>t . (4.8) 

This implies the estimate 



J Wi(t, x, y Q ) dx < 2Qi 4 (t) < 2A Q i 4 *" 1 



for t > to, i G {1,4} and any y G [—1, 1]. An entirely similar argument applied to Q12, Q23, ■ ■ ■ 
yields the estimates 

J w l (t,x,y )dx < 2A t 4e -\ J w l (t,x ,y)dy<2A t Ae - 1 . (4.9) 

for i = 1,2, 3, 4, xq, y$ G [—1, 1] and t > t . 

STEP 2: Pointwise estimates. Using the integral bounds (4.9), we now seek a uniform bound of 
the form 

Wi(t,x,y)<C (4.10) 
for some constant Cq and all x,y G [—1, 1], t > 0. 

To prove (4.10), let t 1— ► x(i), i 1— >• y(i)) G [—1,1] be solutions of 

x = x + 1, y = y + 1 . 

Call 

A(t) = / (toi +t04)(t, x,y(t)) dx. 

Jx(t) 

From our previous estimates (4.9) it trivially follows 

A(t) < AAot 40 - 1 . (4.11) 
The time derivative of A(t) is computed as 
dA I' 1 

— =- (x(t) + l)( Wl +w 4 ){t,x(t),y(t)) + / [OtWi + (y(t) + l)d yWl + d t w A + (y(t) + l)d y w 4 ] dx 

at Jx(t) 

= - (x(t) + l)(wi +W4)(t,x(t),y(t)) - / [wi + (x + l)9 x ^i + ^4 + (x - l)d x w 4 ] dx 

Jx(t) 

= - (x(t) + l)(toi +W4)(t, x(t),y(t)) - / (w\ + Wi)dx 

Jx(t) 

- 2w! (t, l,y(t)) + (x(t) + l) Wl (t, x(t),y(t)) + / Wl dx 

Jx(t) 

-2w 4 (t,l,y(t)) + (x(t)-l)w 4 (t,x(t),y(t))+ I w 4 dx 

Jx(t) 

<[x(t)-l-(x(t) + l)]w A (t,x(t),y(t)) = -2w A (t,x(t),y(t)) . 

9 



This implies 



w 4 (t,x(t),y(t)) < 



IdA 
2~dt 



The total derivative of w\ along a characteristic line is now given by 

(t, x(t),y(t)) =w 2 w 4 - wiw 3 -wi< w 2 w A - w 1 < ^w 2 

k(t) f-dA" 

In turn, for t > t this yields the inequality 
Wl (t,x(t),y(t)) < e -(*- tn ) 



-dA 

IT 



w 1 (t )+ [ e s - tn fc( S )(-^( S ))ds 



<e- ( *- to) [u; 1 (to)+^(to)fe(to)] +e- (t -* o) / A(s)(e s - to k(s))' ds. 

Jto 



(4.12) 



(4.13) 



The first term on the right hand side of (4.13) approaches zero exponentially fast. Concerning the 
second, we have 

e -(*-*°) J* A(s) e s "* n (k(s) + k'(s)) ds < J* e~ {t -^ 2A^ e ~ 1 (d\ns + -^j ds . 

This also approaches zero as t — > oo. Repeating the same computations for all components, we 
conclude that for some time t\ sufficiently large there holds 



Wi(h,x,y) < 



for all x,y G [—1,1]. 



(4.14) 



By continuity, the inequalities in (4.14) remain valid for all x,y in a slightly larger square, say 
[— 1 — e, 1 + e]. For t > t\ we now define 

M(t) = max{wi(t,x,y) ; i = 1,2,3,4, x, y € [-1 - e t_tl e , 1 + e i_tl e]| . 

From the equations (3.2) and (4.14) it now follows 



±M(t)<-M(t) + M 2 (t)<^, 



M(h)<-. 



M(t) < 



< e~ T • e 



for all r > t\. 



Returning to the original variables Ui = e T Wi, this yields 

m < e u 

in a whole neighborhood of the point P* = (t*,x*). Hence P* is not a blow-up point. 
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□ 



5 - A tentative blow-up scenario 

For a solution of the rescaled equation (3.1), the total mass 



f 1 f 1 4 

m(t) = J J ^2wi(t,x,y)dxdy 



may well become unbounded as t — > oo. On the other hand, the one-dimensional integrals along 
horizontal or vertical segments decrease monotonically. Namely, if t <— > y(t) satisfies y = y + 1, 
then 



d * 



dt 

Similarly, if x = x — 1, then 



[wi(t, x, y(t)) + w 4 (t, x, y{t))] dx < . 



d f 1 

— J [w 3 (t,x(t),y) +w 4 (t,x(t),y)]dy <0. 

Analogous estimates hold for the sums w\ + W2 and u>2 + W3. Therefore, a bound on the initial 
data 

w i(0, x, y) E [0, M] for all x,y€[-l,l], 

yields uniform integral bounds on the line integrals of all components: 

J Wi(t, x, y) dx < 4M , J Wi(t,x,y) dy < 4M . (5.1) 



If finite time blow-up is to occur, the mass which is initially distributed along each horizontal 
or vertical segment must concentrate itself within a very small region, thus forming a narrow packet 
of particles with increasingly high density. A possible scenario is illustrated in fig. 3. A packet of 
1-particles is initially located at Pi. In order to contribute to blow-up, this packet must remain 
within the unit square Q. At P2 these 1-particles interact with 3-particles and produce a packet of 
4-particles. In turn, at P3 these interact with 2-particles and produce again a packet of 1-particles. 
After repeated interactions, the packet of alternatively 1- and 4-particles eventually enters within 
the smaller square Q' . After this time, it interacts with a packet of 2-particles at P5 (transforming 
it into a packet of 1-particles) and eventually exits from the domain Q. 

To help intuition, it is convenient to describe a packet as being "young" until it enters the 
smaller square Q', and "old" afterwards. To maintain a young packet inside Q, one needs the 
presence of old packets interacting with it near the points P2, P3, P4 ■ ■ ■ On the other hand, after 
it enters Q', our packet can in turn be used to hit another young packet, say at P5, and preventing 
it from leaving the domain Q. 

As t —>■ 00, the density of the packets must approach infinity. One thus expects that most of 
the mass will be concentrated along a finite number of one-dimensional curves. Say, the packet 
of alternatively 1- and 4-particles should be located along a moving curve 714 (i, #), where 9 is 
a parameter along the curve. The time evolution of such a curve is of course governed by the 
equations 

d d 
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depending on whether 714 (t,9) consists of 1- or 4-particles. The presence of interactions impose 
highly nonlinear constraints on these curves. For example, the interaction occurring in P5 at time 
t implies the crossing of the two curves 714 and 712, namely 

TmM) =7i2M) = P 5 

for some parameter values 9, 9. The complicated geometry of these curves resulting from the above 
constraints has not been analyzed. 
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